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Abstract. In survival or reliability studies, the mean residual life or life expectancy 
is an important characteristic of the model. Here, we study the limiting behaviour of the 
mean residual life, and derive an asymptotic expansion which can be used to obtain a 
good approximation for large values of the time variable. The asymptotic expansion is 
valid for a quite general class of failure rate distributions — perhaps the largest class that 
can be expected given that the terms depend only on the failure rate and its derivatives. 
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the mean residual life 
1. Introduction and Background 



In life testing situations, the expected additional lifetime given that a component has 
survived until time t is a function of t, called the mean residual life. More specifically, if 
the random variable X represents the life of a component, then the mean residual life is 
given by m{t) = E{X — t\X > t). It is well known that the mean residual life is related 
to the survival (reliability) function F by 

1 r°° — 

(1.1) rn{t) = =^J^ F{u)du, 

and to the failure rate (hazard function) r = —F '/ F by 

/OO /' PU ~\ /"OO 

exp I - r{x) dxjdu = e^^^ ^ e'^^") du, 

where 

(1.3) Rit)= [ r{x)dx = - log F{t) 

Jo 

is the integrated failure rate (cumulative hazard function). We also have 

(1.4) m\t) = r{t)m{t) - 1, 

(see Calabria and Pulcini (1987), for example.) 

The mean residual life has been employed in life length studies by various authors, 
e.g. Hollander and Proschan (1975), Bryson and Siddiqui (1969), and Muth (1977). Muth 
(1977) observed that the failure rate takes only the instantaneous present into account, 
whereas the mean residual life takes the complete future into account. Meilijson (1972) 
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has studied certain limiting properties of the mean residual life. A smooth estimator of 
the mean residual life is given by Chaubey and Sen (1999). 

In this paper, we undertake a detailed study of the limiting behaviour of the mean 
residual life O, and derive an asymptotic expansion (Q which can be used to obtain 
good approximations for large values of the time variable. The asymptotic expansion is 
valid for a quite general class of failure rate distributions — perhaps the largest class that 
can be expected given that the terms depend only on the failure rate and its derivatives. 

For the family of age smooth distributions, Rojo (1996) has established, for large values 
of the time variable, a relationship between the mean residual life and the failure rate in 
terms of the index p of regular variation. Calabria and Pulcini (1987) noted a relationship 
between the limiting behaviour of the mean residual life and the failure rate. This rela- 
tionship was developed further by Chaubey and Sen (1999) for the class of distributions 
having nondecreasing failure rate. Our approach provides a considerable improvement on 
their approximation, and moreover does not require that the failure rate be nondecreasing. 



By applying L'Hopital's rule to ()1.1|) . Calabria and Pulcini (1987) derived the relation- 
ship 



2. Limiting Behaviour 



(2.1) 
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provided the latter limit exists and is finite. They then used ()1.4|) to conclude that 
limt^oo "^'(^) = 0, or equivalently, that 



Unfortunately, one cannot infer ()2.2|) from ()2.1|) unless one also assumes that limf^oo r{t) 
is finite and strictly positive. For a counterexample, fix positive constants a and b and 
consider the linear mean residual life m{t) = a + bt with corresponding failure rate r{t) = 
{l + b)/{a + bt). The class of distributions with linear mean residual life has been studied 
by Hall and Wellner (1981) and Oakes and Dasu (1990). Counterexamples that satisfy 
limt^oo ''^(^) = oo also exist, but as these tend to be somewhat more complicated, further 
discussion is deferred to the end of this section. A more detailed study of the limiting 
behaviour of the mean residual life follows. 

Theorem 1. Let r denote the failure rate, and let m denote the corresponding mean 
residual life. For t > define j3{t) = ess sup{r(x) : x > t} and a{t) = ess inf{r(x) : x > 
t} . Then 



(2.2) 



lim r{t)m{t) = 1. 




< liminf m(t) < limsupm(t) < lim 

t^OO t^OO t^OCl 




Proof. Since a (respectively, f3) is clearly nondecreasing (nonincreasing) , 




implies 
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□ 



Corollary 1. The failure rate r and the mean residual life m have limiting behaviour 
related by 

1 1 

-— < liminf m(t) < limsupm(t) < . 

limsupr(t) t^oo i^oo liminir(t) 

Note that Corollary ^ implies ()2.1|) without the assumption that limt^oo^/r{t) exists; 
for example, if limsup^^^ r(t) = 0, then limt^oo^^l^) = oo. Corollary ^ also implies that 
if lim^^oo ''^(^) exists (finite) and is strictly positive, then ()2.2|) holds. 

The limiting reciprocal relationship ()2.2p between the failure rate and the mean residual 
life may be interpreted as an approximation or asymptotic formula m ~ s, where s = 1/r. 
By imposing suitable conditions on s and its derivatives, it is possible to refine this 
approximation by introducing additional terms into the asymptotic formula. We shall 
carry out this programme in ^ Even in cases where the reciprocal relationship ()2.2j) fails, 
one can sometimes obtain reasonably precise information about the limiting behaviour 
of the product of the failure rate and the mean residual life by studying the limiting 
behaviour of s and its derivatives. By putting E = [0, oo) in Corollary |21 below, it can be 
seen that for continuously differentiable s, if limj^oo s'{t) = 0, then ()2.2|) holds. However, 
it may happen that limt^oo s'{t) does not exist. It turns out that if \s'{t) \ is not ultimately 
"too big too often," one can still say a good deal about the product r{t)m{t) when t is 
large: see Theorem0]and Corollary|21below. There are also failure rates for which nonzero 
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values of linit^oo ^'(t) are possible. In such cases, the reciprocal relationship ()2.2|) may 
fail. Nevertheless, we have the following result. 

Theorem 2. Denote the failure rate and its reciprocal by r and s = 1/r, respectively, 
and let m denote the corresponding mean residual life. Let R denote the integrated failure 
rate p.3p . Suppose that 



(2.3) lim s(t)exp(-i?(t)) = 0, 

»oo 



and that 



(2.4) lim(l - s'(t))-^ exists (finite). 



Then 



lim r{t)m{t) = lim (1 — s'{t)) ^ 



t— >oo t— >oo 



Proof. Apply L'Hopital's rule to 

r(t)m(t) = / e-^(") dx / (s(t)e-^W). 



□ 

When applying Theorem |21 in practice, it may often be easier to verify conditions that 
imply the hypothesis ()2.3j) . as opposed to verifying ()2.3|1 directly. We give an example of 
such a condition below. 
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Theorem 3. Denote the reciprocal of the failure rate by s, and let R denote the integrated 
failure rate ()1.3|) . Suppose that lim^^ooess sup{s'(a;) : x > t} < 1 and that s{t) is finite 
for all sufficiently large values oft. Then ()2.Hj) holds. 

Proof. There exist numbers L < 1 and to > such that for all t > to, s{t) is finite and 
ess sup{s'(x) : X > t} < L. For t > to, we have 



r /■* s'(x) 

s{t) exp{-R{t)) = exp <^ logs(to) + / dx 

I J to ^{^) 



t 

— I r(x) dx 

to 



= s(to) exp < / — l)r(a;) dx — / r(x) dx 

I J to Jo 

< s{to) exp < ( ess sup s'{x) — l) / r{x) dx — R{tQ 

< sito) exp {(L - l){R{t) - R{to)) - i?(to)} 
= s(to)e-^^(*«)(F(t))'~''. 



Since limj__oo F{t) = and L < 1, (j23I) follows. □ 

A condition that implies the hypothesis limj^ooGSS sup{s'(x) : x > t} < 1 of TheoremEl 
and is typically even easier to verify is limsup^^go < 1. Thus, Theorems |21 and El 
can be readily applied to a wide variety of situations in which the long term behaviour 
of the rate of change of the reciprocal of the failure rate is known. A very simple yet 
illustrative example is provided by the distribution with fractional linear failure rate 
r{t) = l/(c + dt). Here s' has constant value d, and thus if < 1, then Theorems |21 
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and El imply that limt^ooi^{t)m{t) = 1/(1 — rf). Of course, in this case it is easy to verify 
this fact directly, as the mean residual life function is linear: m{t) = (c + dt)/{l — d). 
However, Theorems El and El are equally applicable to cases in which it may be difficult or 
impossible to determine the long term behaviour of the mean residual life directly. 

The issue of practical usefulness apart, there are compelling theoretical reasons which 
point up the significance of Theorem |21 as well. Recall the observation that ()2.2|) may be 
interpreted as an approximation or asymptotic formula m ~ s. Viewing the conclusion of 
Theorem |21 in the same light yields 

m ~ = s{l + s' + {s'f + {s'f + ■■■), \s'\ < 1, 

1 — s' 

which agrees with the first three terms of the asymptotic expansion ()3.6p . Thus, in 
some sense. Theorem |21 represents a strengthening of Calabria and Pulcini's first order 
result ()2.2j) — subject to the appropriate conditions — to third order. 

We next address the problem of what happens when s' exists but lim^^oo s'{t) does not. 
As might be expected, in general some information about the limiting behaviour of the 
product r(t)m(t) is lost. Nevertheless, in many cases one can at least bound rm in terms 
of the essential supremum of s' . We'll see as a result that even if limsup^..^^ > 0, 

as long as \s'{t)\ is not ultimately "too big too often" then ()2.2p holds. 

Theorem 4. Let r denote the failure rate, and let m denote the corresponding mean 
residual life. Let s = 1/r, and put X{t) := ess sup{|s'(x)| : x > t}. Suppose that r{t) 
is positive and continuously differentiable for all sufficiently large values of t, and that 
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limt^oo (^) < 1- Then for all sufficiently large values oft, 1/(1 + A(t)) < r{t)m{t) < 
1/(1 -A(t)). 

Proof. By Theorem IHl condition ()2.3|) holds. Thus, the hypotheses permit us to integrate 
by parts and discard the hmit at the upper range of integration in the integrated term. 
We have for all sufficiently large values of t, 



r(t)m(t) = r(t)e^(*) / r{x)e~^^''h{x) dx = I + r{t)e^^'^ e-^^'^h' {x) dx. 



But 



e-^^'^h'ix) dx 



< ess sup \s'{x)\ / e'^^""^ dx. 

x>t Jt 



Thus, for all sufficiently large values of t. 



1 - \{t)r{t)m{t) < r{t)m{t) < 1 + \{t)r{t)m{t). 



□ 

It follows that the reciprocal relationship ()2.2|1 may hold even if limsup^^^ > 0, 

as long as lim„^oo 'S'(tri) = for "most" sequences ti < t2 < ■ • ■ — ^ oo. Here, "most" is in 
the sense of Lebesgue measure. The following result makes this observation more precise. 



Corollary 2. Lets denote the reciprocal of the failure rate. Suppose that s{t) is finite and 
continuously dijjerentiable for all sufficiently large values oft. Suppose further that there 
exists a subset E of the interval [0, oo) whose complement in [0, oo) is of Lebesgue measure 
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zero, and such that for every sequence ti, ^2? ■ ■ ■ of elements of E with lim„^oo tn = oo we 
have lim^^oo ■s'(tn) = 0. Then ()2.2|) holds. 

Proof. As customary, denote the indicator function of a set A by xa- Suppose that 
E C [0, oo) satisfies the hypotheses of the Corollary. Then the complement of E in [0, oo) 
has Lebesgue measure zero, and Imix^oo s' {x)xe{.x) = 0. Therefore, if e > is given, 
there exists a suitably large value of t such that \s'{x)x[t,oo)(yX)\ < e for all x E E. By 
definition of the essential supremum, this implies (in the notation of Theorem 0)) that 
\{t) = ess sup{|s'(x)| : X > t} < e. As A is a nonincreasing function and e > is 
arbitrary, it follows that lim^^oo A(t) = 0. □ 

We conclude this section with an example of a distribution in which limf^oo^(^) = oo, 
but limt^oo '''('^)"^(''^) 7^ 1- Let a,b,c,d be positive constants satisfying a > b, d > 2b, and 
> (a + b)d. Consider the mean residual life defined by 



m{t) 



a + 6sin(t^) 



t > 0. 



c + dt 



We have 



m'{t) 



2btcofi{f) {a + bsi\i{t^))d 
c + dt (c + dt)2 



Thus, 



lim inf m'{t) 



2b/d, 



limsupm'(t) = 2b/ d, 



and hence 



lim inf r(t)m(t) = 1 + liminfm'(t) = 1 — 2b/ d < 1, 
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whereas 



lim sup r(t)m(t) = 1 + limsupm'(t) = 1 + 2h/d > 1. 

t—KX> t— »00 



Since 



{c + dtf{l + m'{t)) = {c + dtf + {c + dt)2btcos{f) - {a + bsm{f))d 
>c^ + 2cdt + dH"^ - 2bct - 2hdt^ -{a + b)d 
= (d - 2b)dt'^ + {d- b)2ct + - (a + b)d 

is clearly positive for t > 0, it follows that r{t) = (1 + m'{t))/m{t) > for i > 0. 
Furthermore, 

c+ {d + 2bcos{f))t d 
^ ~ a + 6sin(t2) c + dt 

implies \imf_^^r{t) — oo. 

3. Asymptotic Expansion 

Under certain conditions, the mean residual life has an asymptotic expansion in terms 
of the failure rate and its derivatives. An initial attempt in this direction was made 
by Chaubey and Sen (1999) for the class of distributions having nondecreasing failure 
rate. However, it is easy to envision situations where, say with regular maintenance, even 
an ultimately nondecreasing failure rate may be an inappropriate model. Therefore, we 
provide an alternative approach that requires no monotonicity assumptions on the failure 
rate. 
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We take Chaubey and Sen's asymptotic formula 

1 r'{t) r"(t) \ 

m(t) = -rr - , + Q L , t ^ oo, 

(with error term corrected), as a point of departure. Unfortunately, their derivation is not 
rigorous, and as a result, certain growth conditions on r and its derivatives are omitted. 
Nevertheless, their approach can, in principle, be made to work, and thus one can show 
that under suitable conditions on the failure rate there is an asymptotic series development 
to arbitrary order that begins 

mr^r-^- r'r-'' - r"r-^ + (3(r')' - r'") r'' + (lOrV" - r"") r"^ + • • • , t ^ oo. 

Here, we have made the abbreviations m = m{t), = l/r{t), etc. More explicitly, for 
each positive integer n, we have 

n-l 

(3.1) m ~ ^^Cfcr"*^"^ + (additional terms), 

fe=0 

where as i — > oo, the additional terms tend to zero more rapidly than Ckr~''~^ for 1 < 
k < n — 1. Here, Cq = 1, Ci = 0, C2 = — r', C3 = — r", C4 = 3(r')^ — r'", and in general, 
is a polynomial in r', r", . . . , r^'^"^^ for k >2, the explicit form of which is given by 

(3.2) 5: M.^nj^Girir ■ 

p=0 j>l J ' I / 

where the inner sum is over all nonnegative integers ai, q;2, . . . such that X^j>i OLj — p and 

Ej>i(j + i)«j = 
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Rather than attempt a rigorous proof of ()3.1|) and ()3.2|) along these hnes, we instead 
develop an alternative, operator-theoretic approach that appears to be simpler, yet more 
powerful. By way of motivation, observe that since we are interested in the situation 
when ()2.2j) holds, it makes sense to let s = 1/r and rewrite the differential equation ()1.4j) 
in the form (1 — sD)m = s, where D denotes the derivative operator. Abbreviating sD 
by O, we might hope to find a meaningful way to define the inverse operator (1 — 0)^^ 
in such a way that given suitable growth conditions on s and its derivatives, we have a 
legitimate asymptotic expansion 

m = (1 - Q)-h ~ s + Q{s) + Q\s) + Q^{s) + ■■■ . 

That this is indeed the case is the main content of Theorem El in which an explicit formula 
is also given for Q^{s), the kth term of the expansion. A significant step in this direction 
is provided by the following result. 

Lemma 1. Let s denote the reciprocal of the failure rate, and suppose that for some 
nonnegative integer n, s is n + 1 times continuously differentiahle on the positive real line. 
As usual, let R denote the integrated failure rate ()1.3|) . and let D denote the derivative 
operator. For positive integers k, let := (sD)^ and for convenience, let G° := s. If 
in addition, we have lima;^oo exp(— = for each k = 0,1, ... ,n, then the 
mean residual life can be expressed in the form 

^ POO 

(3.3) m{t) = sit) + J2 + e"^^*^ / r(x)e-^(")e"+^(s(x)) dx, t > 0. 
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Proof. The conditions on s permit us to integrate by parts n + 1 times and discard the 
Umits at the upper range of integration. Thus, integrating by parts once, we have 



(t) = e^W / r(x)e-^(")s(x) dx = s{t) + e^W / r(x)e"^(")0(s(x)) dx. 



It is now evident that the claimed formula ()3.3|) can be proved using mathematical induc- 
tion and repeatedly integrating by parts. 

Alternatively, note that (1 — s{t)D) exp(i?(t)) r(x) exp{—R{x))h{x) dx = h{t) for all 
locally integrable functions h for which the integral converges, and so we may define 

POO 

(1 - Q)-^h{t) = e^W / r(x)e-^(")/i(a;) dx 



for all such functions h. If we now take h = s and then h = G""'"^s, then in light of the 
identity 

(1 - Q)-^ = 1 + e + 02 + . . . + 0" + (1 _ 

we have 

/oo 
r{x)e-^^''ls{x)dx 

n 
k=l 

^ POO 

= sit) + ^^(sit)) + e^(*) / r(x)e-^(")e"+^(s(x)) dx, 
as claimed. □ 
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Theorem 5. Let s denote the reciprocal of the failure rate. Suppose that for some non- 
negative integer n, s is n + 1 times continuously difjerentiahle on the positive real line, 
and that liiiit^oo ■s(t) exp(— /^(t)) = 0, where as usual R denotes the integrated failure 
rate (jl.Hj) . Let D and 6 he as in LemmaUi Suppose that for each k = 0,1, . . . ,n we have 
s(t)s(^+^)(t) = o(sW(t)) as t ^ oo. Then e'=+^(s) = o(0'=(s)) for all < k < n, and we 
have the asymptotic expansion 

n 

m{t) = + ^e*=(s(t)) + o(e"(s(t))), t^oo 

k=l 

for the mean residual life. The terms Q''{s) may be expressed in terms of s and its 
derivatives by means of 

k 

(3.4) e'=(.) = . Yl d{j„...,j,)l[s'~^^\ 

jl,--;jk>0 P=l 

where the sum extends over all nonnegative integers jp [1 < p < k) and the coefficients 
d{ji, . . . ,jk) are generated by the polynomial equation 

k k p 

(3.5) Yl d{n,...,j,)l[xir = m2''^ 

ji,-,jk>0 p=l p=l j=i 

in the indeterminates Xi, . . . ,Xk. 

Proof. The stated formula for 0*^(5) is a special case of a more general result. See 
Snyder (1982). We next show that, under the hypotheses on s, Q^^^{s) = o(G^(s))) for 
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< k < n. Since O'^+Hs) = sDQ''(s), we have 



k 

jl,--;jk>0 P=l 



k 

s'Q\s) + s Yl d{n,...,jk)sDl[{D^-s) 

iivjfe>o p=i 

k k 



Note that the generating function ()3.5|) imphes that ji + ■ ■ ■ + jk = k above. In particular, 
since < k < n, each jq < n, and so each sD^'J^^s = o{D^''s). Also, s' = o(l) because 
ss' = sDs = o{s). Therefore, we have 

/ k k \ 

^(7=1 7l,...,7t>0 P=l ^ 



9=1 jl,-;jk>0 

o{e\s))+o{ke\s)) 



o 



Next, observe that since limj^oo ■s(t) exp(— _R(t)) = and 6'^+^(s) = o(6'^(s)) for all 
/c = 0, 1, . . . , n, we have limj^oo exp(— i?(t))6^(s(t)) = for all A; = 0, 1, ... , n. Therefore, 
in light of Lemma ^ it remains only to show that 



" POO 

m(t)-^0'=(s(t)) = e^W / r(x)e-^(")e"+i(s(x))rfx = o(e"(s(t))), t ^ oo. 

k=0 
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Since = 0(0" (s)), Lemma ^ gives 

m(t) - J2 ©"'(sW) = / r(x)e-^(")e"+^(s(x)) dx 
k=o 



o(m(t)-Xi0'(^W))- 

^ fc=0 ^ 



Therefore, 

n-l 



(1 - 0(1)) (m(t) - J2 = 



from which it follows that 



n-l 
k=0 



or equivalently, 



n 

m{t)-Y,Q\s{t)) = o{Q%sm 

k=0 



as claimed. □ 



/// 



We have e%s) = s, 0(s) = ss', e\s) = s{s'y + s^s", Q\s) = s{s'Y + 4sh's" + sh 
and e\s) = s{s'y + lls^s'^s" + As^{s"f + 7sh's"' + s^s"". Thus, under the conditions 
of Theorem 13 our asymptotic expansion begins 

(3.6) m^s + ss+ s{s')^ + s^s" + s{sf + As'^s's" + sh'" + ■ ■ ■ , t ^ 00. 
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